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Abstract 

A "persistence" exponent 6 has been extensively used to describe the nonequilib- 
rium dynamics of spin systems following a deep quench: for zero-temperature homo- 
geneous Ising models on the d-dimensional cubic lattice Z'^, the fraction p(t) of spins 
not flipped by time t decays to zero like t~^^'^^ for low d; for high d, p{t) may decay 
to p(oo) > 0, because of "blocking" (but perhaps still like a power). What are the 
effects of disorder or changes of lattice? We show that these can quite generally lead 
to blocking (and convergence to a metastable configuration) even for low d, and then 
present two examples — one disordered and one homogeneous — where p{t) decays 
exponentially to p(oo). 



In modelling the nonequilibrium dynamics of spin systems following a deep quench, the 
following question naturally arises [jl], |], ^ ||, ^ given a spin system at zero tem- 
perature with random starting configuration and evolving according to the usual Glauber 
dynamics, what is the probability p{t) at time t that a spin has not yet flipped? 

For the homogeneous ferromagnetic Ising model on Z'^, this probability has been found 
to decay at large times as a power law p{t) ~ t-^^'^) m, |, I for < 4. The "persistence" 
exponent 6{d) is considered to be a new universal exponent governing nonequilibrium dy- 
namics following a deep quench ||^. The persistence problem can be extended to positive 
temperatures by considering the dynamics of the local order parameter rather than that of 



single spins M, OO, 11 



In this paper we confine our attention to dynamics at zero temperature in infinite spin 
systems. In the usual case of asynchronous updating, a spin is chosen at random (this 
can be made precise for infinite systems, as in |1^) and then: always flips if the resulting 
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configuration has lower energy, never flips if the energy is raised, and flips with probability 
1/2 if the resulting energy change is zero. We will consider these dynamics for random 
initial configurations (in which each spin is equally likely to be up or down, independent 
of the others) in both disordered ferromagnets and spin glasses with continuous coupling 
distributions, and also for uniform ferromagnets on lattices other than Z'^ (e.g., hexagonal 
lattices in 2D). 

Our first result is that the persistence phenomenon as described above is unstable to 
the introduction of randomness into the spin couplings, or even to some changes in lattice 
structure. For the random ferromagnet, spin glass, 2D hexagonal ferromagnet, and others 
to be discussed below, we will see that a positive fraction of spins never flip and every spin 
flips only finitely many times. 

The "frozenness" of a nonvanishing fraction of spins (sometimes referred to as "blocked" 
spins 0) has been reported in numerical simulations of Ising ferromagnets on Z'^ with d > 4 

and g-state Potts models on square lattices for g > 4 ||^. The problem can then be 
recast by restricting attention to only those spins that eventually do flip, and asking for the 
conditional probability that such spins haven't yet flipped by time t. Simulations of Potts 
models appeared to indicate that this probability (proportional to p(t) — p{oo)) also decays 
as a power law at long times [§] (however, some curvature on their log- log plots was noted). 

We will examine the same question for disordered Ising systems and also for homogeneous 
systems that show blocking. Although we cannot yet answer this question in general, we will 
present calculations on two systems, the homogeneous ferromagnet on a quasi- ID "ladder" 



and the ID disordered spin chain |T^, showing that p{t) — p{oo) decays exponentially as 
t oo. Exponential decay for d > 2 will also be discussed. 

Persistence and local nonequilihration. The analysis of persistence exponents suggests 
that the fraction of sites that remain in the same phase (for T > 0) or spin value (at T = 0) 
from time ti to time ^2 tends to zero for 1 <C ti -C ^2- It therefore implies the presence of 
local nonequilihration (LNE) that in any fixed, finite region, there exists no finite time 
after which the spins within remain in a single phase; that is, domain walls forever sweep 
across the region. At zero temperature, the presence of LNE means every spin flips infinitely 
often (in almost every sample). 

Why does the decay to zero of p{t) (coming from the analysis of persistence exponents at 
T = 0) suggest that every spin flips infinitely many times? Suppose instead that a positive 
fraction of spins flip only finitely many times. Then it is reasonable to expect that a (smaller 
but still positive) fraction of spins never flip, and p{t) would not decay to zero. While not 
proved in general, this argument applies to all systems treated here. 

It was proved in 0] (see also |jl2|) that, in the homogeneous Ising ferromagnet on the 
square lattice with a random initial spin configuration, every spin indeed flips infinitely 
often at zero temperature, consistent with persistence results in the literature. (Similar 
results apply to several other systems, and can be extended to positive temperature with 



the local order parameter in a region replacing individual spins |12]. 



Blocking. What about the zero-temperature dynamics of systems with continuous dis- 



2 



order? In any dimension and on any lattice, it can be proved for these (and many other) 
systems that every spin flips only finitely many times. These systems exhibit "blocking" and 
for them p{t) does not decay to zero. 

These are examples of a general result [|14| applying to the dynamical evolution (following 
a deep quench) of infinite- volume Ising spin systems with Hamiltonian 

= — X! Jxy(^x<^y , (1) 

<xy> 

where the sum is over nearest neighbors. If the distribution of couplings is continuous with 
finite mean, then it can be proved that every spin flips only finitely many times (for almost 
every cr°, realization uj of the dynamics and realization J' of the couplings). 

The proof of this theorem yields a more general result that shows that, even without 
the continuity assumption on the distribution of couplings, for almost every J', a^, and u, 
there can be only finitely many flips of any spin that cause a nonzero energy change. This 
is why the above result applies to ordinary spin glasses and random ferromagnets with a 
continuous distribution of couplings (e.g., Gaussian or uniform): the probability of a "tie" 
in any sum or difference of a given spin's nearest-neighbor coupling strengths (and therefore 
the probability of a spin flip costing zero energy) is zero, and the result follows. 



We sketch the proof here; for further details, see [Q. Let a* be the value of at time 
t for fixed u, and J'. Let 



E(t) = -(l/2) J2 Jxyaial (2) 

y:\x-y\=l 

where the bar indicates an average over J', a^, and uj. By translation-ergodicity of the 
distributions from which J', cr°, and cu were chosen, and using the assumption that \Jxy\ < oo, 
it follows that E{t) exists, is independent of x, and equals the energy density (i.e., the average 
energy per site) at time t in almost every realization of JT", and u. 

Because every spin flip lowers the energy, E{t) monotonically decreases in time (note that 
E{0) = 0) and has a finite limit E{oo) (> ~d\J.j.y\). Now choose any fixed number e > 0, 
and let be the number of spin flips (over all time) of the spin at x that lower the energy 
by an amount e or greater. Then — oo < E{oo) < —eN^ so that for every x and e > 0, A''^ is 
finite. Let be the minimum energy (magnitude) change resulting from a flip of cTx] then 
although ex varies (differently in each J) with x, it is sufficient that it is strictly positive. 

This result applies also to homogeneous systems on certain lattices, such as Ising ferro- 
magnets on lattices with an odd number of nearest-neighbors so that ties in energy cannot 
occur. Such lattices include the hexagonal (or honeycomb) lattice in 2D, and the double- 
layered cubic lattices Z'^x {0, 1} (i.e., a "ladder" when d = 1, two horizontal planes separated 
by unit vertical distance when d = 2, and so on) |]15| . 



As for blocking in these systems, it is elementary to show that a positive fraction of 
spins will never flip. Consider first the hexagonal lattice. If the spins on any single hexagon 
are all up or all down, they will form a stable configuration that will never change. Such 
configurations (and of course similar larger-scale ones) occur with positive density in almost 
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every cr°. Similarly, in the ladder, any square with all spins up or down is stable. The 
extension to general Z"^ x {0, 1} is straightforward. 

Turning to disordered systems, consider first the random ferromagnet on . For almost 
every JT", there will be a positive density of plaquettes whose couphngs satisfy the following: 
on each of the four corners, the sum of the two couplings that connect to adjacent corners 
of the square is greater than the sum of the two couplings to sites outside the square. If the 
spins at the four corners initially are all up or all down, the spin configuration on the square 
will again be stable. A similar construction can be used for general d and for spin glasses. 

To summarize, our first result has been to prove that many ordered and disordered spin 
systems display two important zero-temperature dynamical properties, which when taken 
together lead them to exhibit novel persistence behavior. The first concerns the presence 
of blocking, meaning a positive fraction of spins never flip. In the systems we treat, this 
is a zero time property in that some of the spins are blocked by the nature of cr^ (and 
J\ regardless of the dynamics realization. The second property concerns infinite time: the 
existence of a limiting (metastable) spin configuration 0"°°, since every spin flips only finitely 
many times. Although the second of these properties probably implies the first, the first 



does not imply the second |jT6[. Our next result is to show that for at least some of these 
systems, these two properties lead to an exponential (as opposed to power law) decay of the 
quantity p(t) — "pipo) at large times. 

Exponential decay. In this section we study the large-time behavior of p(t) — p(c)o), the 
probability that a spin will flip at some time but has not yet flipped by time t. We will 
prove that this quantity decays exponentially by showing the same for the larger probability 
p(t) that a spin will flip at some time after t (whether or not it has flipped before). We 
consider two systems, one homogeneous (the uniform Ising ferromagnet on the ladder) and 
one disordered (the \D continuously disordered Ising chain). 

Consider first a homogeneous system where every site has an odd number M of neighbors. 
(Systems such as ±J spin glasses where the signs are disordered but not the | J^j/I's also fall 
into this category of examples.) Consider at time r, all sites y such that the spin at y will flip 
after time r, and denote by Cxij^ the cluster of such sites that contains x (an empty cluster 
if the spin at x will not flip after time r). We will show below for the ladder model that with 
r = 0, the distribution of the number of sites |Cj;(r)| in these clusters has an exponential 
tail; i.e., the probabilities for large cluster sizes are bounded by: 

P(|C.(r)| >n) < Ae-^" (3) 

for some A < oo and A; > 0. We next show that this implies exponential decay of p(t). 

Since each flip in Cxij) lowers the energy of that cluster by at least 2 and since the 
total energy of the cluster lies somewhere between —M\Cx[j')\ and M|Ca;(r)| (we take J = 1 
here), it follows that the entire cluster must reach its final configuration after no more than 
^|Cx(t)| flips. Let Ti denote the (random) amount of time after r until the first flip in 
C^.(r), T2 the amount of time after r + Ti until the second flip, etc. Clearly, as long as flips 
are possible, the Tj's are bounded above by independent exponential (mean one) random 
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variables T/. Thus the time of the last flip of x is bounded above by r + T{ + ■ ■ ■ + T'j^^^^f^^-^ 
and so for t > r. 



Pit) - p{oo) < Pit) < Y: PilC.ir) I = n)P{T[ + • • . + TX,„ > t - r) . (4) 

n=l 

The probability density of T{ + ■ ■ ■ + Tj is /(s) = s^~^e'~^ / {j — 1)! and so for t > t, 

^ POO 

p{t) - p{oo) < Y ^e"^" / [s^'^''-^/{Mn - l)\]e-'ds 

n=l -^^-^ 

POO °^ 

< / ^A(e-'=/*Ois'"V(j-l)!]e-^(is 
= Ae-'^/'' r expie-'^^'s - s)ds 



lt~T 

\l„-k't 



= A'e-''' (5) 

where the constants A' and k' depend on A, k, M and r. 

It remains to show that (^) is valid in the ladder ferromagnet with r = 0. (Similar 
arguments work for the ladder antiferromagnet or ± J spin glass.) To do this we note that a 
single plaquette has an initial probability po = 1/8 of having its four corner spins all plus or 
all minus (we call such a blocked plaquette "frigid"). The lattice's sites take integer values 
{xi,yi), with —oo < Xi < oo and yi = 0, 1. A lower bound for the initial number of frigid 
plaquettes can be obtained by considering only those plaquettes whose left edges occur at 
even Xi (we deflne the location of a plaquette by the position of its left edge); such plaquettes 
do not overlap and so their probabilities of being frigid are independent. If the plaquette 
at the origin is frigid, and P(0 2n) is the probability that there is no frigid plaquette 
between and X2n, then 

P(0 2n) < (1 - PoT = exp{-kn) , (6) 

where k = \ log(l — po)\. So the ladder is broken up into flnite segments, bounded to either 
side by a frigid plaquette, whose length distribution has an exponential tail. This yields 
Eq. (I). 

Our second example is a disordered ID spin chain in zero fleld. The analysis is essentially 
the same for either the spin glass or the random ferromagnet, so for specificity we study a 
ferromagnet whose couplings Jz = Jz,z+i are independent random variables taken from the 
uniform distribution on [0,1]. The key idea here is that the chain breaks up into finite, 
disjoint "infiuence segments" whose union is the infinite chain. An infiuence segment is a 
dynamical construct defined (for a given J') as follows: two sites x and y belong to the same 
infiuence segment if and only if either the state of can dynamically induce a change in 
the state of ay or vice- versa (or both). To illustrate, suppose that the coupling is larger 
than both Jx-i and Jx+i', i.e., Jx is a local maximum. Then it is clear that the state of ax 
can be dynamically infiuenced by ax+i but not by ax-i (and similarly, the state of ax+i can 
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be dynamically influenced by cTx but not by That is, no state of the spin (Tx-i can 

alter the sign of the energy change A?i^ that would result from a flip of cr^. To summarize, 
two sites X and y are defined to be in the same influence cluster if and only if either ax can 



influence o"y, or ay can influence cr^, or both [|I7 



Influence segments for the disordered ID chain are then constructed as follows [|I4 
Consider the doubly infinite sequences x = Xm of sites where Jx is a local maximum and 
y = Um ^ {xm,Xm+i) where Jy is a local minimum: the couplings are strictly increasing 
from Hm-i to Xm and strictly decreasing from x^ to i/m- The set of spins at the sites 
{Um-i + 1, Vm-i + 2, . . . ym} determines a single influence segment. 

To see this, note that the spin at ym cannot influence the one at ym + ^ (or vice- versa); 
similarly, at the other end ym-i + 1 cannot influence ym-i- Now consider the spins at Xm 
and Xm + 1, which are within the interval {ym-i + 1, Vm-i + 2, . . . ym}- Clearly, the spin at 
Xm — 1 can never influence the spin at Xm, and the spin at Xm + 2 can never influence the one 
at Xm + 1. So once the spins at Xm and Xm + 1 agree (either initially in 0"° or after either spin 
flips) the final value of every spin in the interval {ym-i + 1, Vm-i + 2, . . . ym} is determined 
through a "cascade" of influence to either side of {xm,Xm + 1} (which is put into effect as 
the Poisson clocks successively ring) until ym-i + 1 and ym, respectively, are reached. 

Given this, the analysis leading to Eqs. (P and (^ applies as before. One needs only an 
estimate analogous to Eq. (|^) for the probability distribution of influence segment sizes. In 
fact, the decay here for large size n is faster than in Eq. (^; the probability of n independent 
coupling random variables being ordered so as to have a single local maximum falls off as 
l/n\ (times an exponential factor). 

In these two examples different factors determine the distribution of dynamical cluster 
sizes: for the homogeneous ferromagnet on the ladder, they're determined by the initial spin 
configuration cr°; for the disordered ID chain, they're determined by the coupling realization 
J. 

In this section we considered two examples, one ordered and one disordered, but both 
one-dimensional (or quasi-one-dimensional). There is another system that shows the same 
behavior in any dimension: the highly disordered spin glass (or ferromagnet) [0, |19|, 
Using similar arguments, this system can also be shown to display an exponential decay 



to its final state [21|. We expect that a related model, in which coupling magnitudes are 



"stretched" in the manner of references |T8|, |T9|, |2y] but only up to a finite length scale, would 
show similar behavior. This last model is of interest because its thermodynamic behavior is 
expected to be similar to that of the ordinary spin glass (or random ferromagnet). 

Discussion. Most work on persistence at zero temperature has examined systems, such as 
the homogeneous Ising ferromagnet on Z'^ in low d, where the quantity p{t) decays to zero as 
a power law. We have shown here that there is a second class of models in which p{oo) > 0: 
these include systems with continuous disorder and homogeneous systems on other lattices. 
In several of these the persistence decay is exponential rather than power law. It would be 
of interest to see whether this fast decay holds in other systems in this general class, such as 
the 2D homogeneous ferromagnet on a hexagonal lattice or an ordinary spin glass with 
d>l. 
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Although we don't know the answer in the general case, we can speculate using a rough 
argument that the answer may be yes. If every spin flips only flnitely often, as time progresses 
an increasing number of spins will "freeze"; i.e., they cease to flip. It is reasonable to expect 
that after some flnite time "unfrozen" spins no longer percolate, so that the dynamics is 
conflned to noninteracting flnite clusters, as in the examples treated here. Of course, there 
remain serious gaps: this is not independent percolation, and the dynamics in the localized 
clusters (should they exist) would need to be worked out, so the conclusion should be treated 
with caution. 

There is a third class of systems not discussed here; in these, a positive fraction of spins 
flip inflnitely often and a positive fraction flip only flnitely many times. One such system is 
the two-dimensional ±J spin glass |T6|. Although it appears that p{oo) > 0, determining 
the large-time behavior of p{t) — p{oo) remains an open problem. 
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